Red solitons: evidence of spatio-temporal instability in spatial-soliton dynamics 
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In x' 2 ' three-wave mixing, the noise-seeded spatio-temporal modulational instability has a dra- 
matic impact on the spatial-soliton dynamics, leading to the counterintuitive observation of a soliton 
with no apparent participation of the high-frequency field in the process. 



Optical spatial solitons are ideally monochromatic 
light beams whose linear diffractive spreading is compen- 
sated by a suitable non-linear phase shift. In x' 2 ' three- 
wave mixing, multi-color solitary beams are sustained by 
energy exchange among the three waves, whose phase- 
velocity and wavefront-curvature mismatch give rise to 
an intensity-dependent phase accumulation 2]. Since 
their first observation in a second-harmonic generation 
(SHG) scheme |3j, X^ spatial solitons have been in- 
vestigated in a number of different experimental condi- 
tions, i.e., in up and down conversion, in bulk, wave- 
guides and periodically-poled crystals, with gaussian and 
vortex-type pumps; their unique steering, addressing and 
switching properties have been studied in great detail for 
potential applications in all-optical interconnects (see ref- 
erence for an updated review on the topic). 

With the exception of two recent experiments, where 
fairly short pulses and long samples were used 0> @j 
all the above mentioned phenomenology has been inves- 
tigated in conditions in which material chromatic dis- 
persion was assumed to play a negligible role. How- 
ever, no matter how long the input pulses are {i.e., how 
narrow the input pulse spectra are), one should expect 
appreciable frequency broadening to occur, due the ef- 
fect of the noise-seeded modulational instability (MI) 
0- In the spatial domain, sizeable MI impact has al- 
ready been shown to cause spatial filamentation of the 
SHG-generated temporal soliton, for operation just above 
threshold ||. By analogy, one would expect temporal 
pulse fragmentation to affect the spatial-soliton regime. 
This process is, however, difficult to detect directly since 
it would require ultra-high temporal resolution in single- 
shot acquisition mode. 

In this paper we present the first evidence for the effect 
of noise-seeded spatio-temporal (ST) MI on x^ spatial- 
soliton dynamics, which we predict to cause a temporal 
breakup on the 10-fs scale. We claim ST MI to be the 
mechanism that explains our counter-intuitive observa- 
tion of a spatial soliton for which self-trapping only oc- 
curs between the low-frequency waves {e.g., the signal 
and idler), with no apparent participation of the high- 



frequency wave {e.g., the pump). This is the reason we 
use the term "red soliton" to indicate such a process even 
though, strictly speaking, no solitons are formed. 

Our first experiment was performed in the regime of 
optical parametric amplification (OPA) of the vacuum- 
state fluctuations, with the usual scheme assumed to lead 
to spatial-soliton formation a pump pulse (527 nm, 
1.3 ps, 32 /im spot FWHM) was launched on the input 
face of a 15mm-long lithium triborate (LBO) crystal, 
operated in non-critical, type-I phase matching. For the 
given pump frequency (f2 p ), we tuned the signal and idler 
frequencies (f2 s and f2j) by changing the crystal temper- 
ature. Since the chromatic dispersion {i.e., the group- 
velocity mismatch GVM and the group-velocity disper- 
sion GVD) experiences large changes with the tuning pa- 
rameter D, = Q s /ilp, the setup allowed easy monitoring 
of the effects of dispersion on the non- linear dynamics. 
We should mention that here that the GVM and GVD 
lengths are respectively 2-4 and 10 2 -10 3 times larger 
than the crystal length, across the whole tuning range. 

Two different diagnostics were realized. The first is 
a traditional, CCD-based, beam-profile detection, which 
allowed recording of the flucnce distribution of the pump 
and signal waves at the crystal output face. Note that 
the idler was not recorded, its photon energy being too 
low for the silicon detector. The second setup provided 
three-dimensional {i.e., ST) mapping of the pump inten- 
sity profile, taken at the same crystal-output plane. The 
scheme is based on an ultrafast non-linear gating tech- 
nique realized by mixing the unknown wave packet 
(WP) with a 200-fs plane pulse generated by nonlinear 
compression of a portion of the same pump-laser pulse. 
The achieved temporal resolution is here of 200 fs. 

Figure ^ shows the fluence profiles along the beam 
diameter of the pump (full line) and of the signal field 
(dashed line) just above threshold for soliton formation. 
Figure^ refers to a crystal temperature T = 115° (which 
corresponds to = 0.684 and X s = 770 nm), close to the 
edge of the tuning curve; Fig.^i refers to T = 140° (and 
thus f2 = 0.620 and X s = 850), i.e. much closer to degen- 
eracy. By looking at the result in Fig. [IJi, one would rea- 
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FIG. 1: Transverse fluence distribution for the output sig- 
nal and pump beams for OPA excitation and two different 
frequency-tuning conditions (a and b), with the same pump- 
pulse energy, E p = 2.1 J. 



sonably affirm that spatial solitons are formed. In fact, 
mutually trapped pump and signal beams of comparable 
peak fluence appear. The pump FWHM beam diame- 
ter (15 fim) is five times smaller than the value obtained 
in the linear propagation regime. In contrast, Fig. ^ 
contains quite a puzzling result. In fact, while the sig- 
nal keeps its narrow size, the focused pump beam al- 
most disappears into the weak background field. The 
background profile perfectly matches the pump-beam 
pedestal in Fig. Note the growth of the signal peak 
fluence from Fig.^ to^>, its value being 20 times larger 
than that of the pump in the second case. By taking 
a crystal of twice the length we obtained very similar 
phenomenology, the diameters of the output self-focused 
beams being only 30% larger than in the previous case. 

A deeper insight into this rather counterintuitive ob- 
servation is given by the results presented in Fig. El (left), 
where the experimental ST intensity distribution of the 
pump WP is plotted for two different f2. The results 
show that: (i) the trailing edge of the pulse (about 200 fs) 
does not experience spatial trapping owing to its low in- 
put intensity. This portion accounts for the pedestal and 
diffracted background in Figs. QJi and respectively. 
Note the slight temporal asymmetry of the entire ST dis- 
tribution, caused by GVM. (ii) The central and leading 
portions, in contrast, always give rise to a spatially self- 
trapped profile. Its peak intensity, however, decreases 
dramatically as the tuning condition approaches degener- 
acy, which explains why the trapping seems to disappear 
in the time-integrated profiles. Time-resolved measure- 
ments closer to degeneracy were prevented by the low 
sensitivity of our CCD video trigger. 

In order to provide a more accurate description of the 
underlying process, we performed high-resolution (4-fs 
step) numerical experiments by solving the three-wave 
parametric equations within the framework of a 3D+1, 
ST model with radial symmetry and temporal dispersion 
up to third order. Figure (center) shows the calculated 
pump profiles for conditions identical to those in Fig. [21 
(left). Here the resolution of the numerical-data presen- 



tation is reduced to 200 fs in order to match the resolu- 
tion of the experimental diagnostic. The excellent agree- 
ment between experiment and calculation proves that the 
key features of the observed phenomenon are determined 
by the \^ process only. Importantly, the full resolu- 
tion plots, presented in Fig. [21 (right), indicate the occur- 
rence of a temporal breakup on the 10 fs scale (the cor- 
responding autocorrelation shows a 25 fs FWHM peak). 
The same structure occurs on the signal and idler waves, 
which appear spatially and temporally self-trapped with 
the pump. As to the dependence on fi, the idler behaves 
in the same way as the signal wave. 

Evidence for the fundamental role played by the di- 
mensionality in the non-linear process is given by the 
plots in Fig. El where results of the experiment, of the 
spatio-temporal, spatial and temporal models are com- 
pared. The results indicate that while the ST model de- 
scribes the experiment fairly well, both the monochro- 
matic and plane-wave approximations completely fail to 
reproduce the large variation of the fluence ratio with the 
tuning parameter, giving an almost constant value (about 
unity) over the entire tuning range. This unequivocally 
proves that the red solitons appear as a genuine ST ef- 
fect. Note that, since the temporal dynamics occurs on a 
time scale much shorter than that of the pulse envelope, 
one should not expect any appreciable effect of the input- 
pulse duration on the observed phenomenology. Direct 
numerical test with 5 ps long pulses did not produce any 
detectable deviation from the results shown in Fig. 

Our claim is that the red solitons appear as a conse- 
quence of the ST MI. In justifying this, we should first 
note that no theory is yet available concerning the in- 
stability of an intense, CW and focused pump beam. 
Therefore we must rely, in first approximation, on the 
predictions valid for the CW and plane- wave regime. In 
this case, quantum-noise parametric amplification has al- 
ready been shown to appear as a branch of the MI pro- 
cess 01 ■ When the ST features of the noise are ac- 
counted for, one should not expect spatial and temporal 
breakup to occur as independent processes. In fact, via 
the phase-matching constraint, the parametric gain (or, 
equivalently, the MI) couples space and time by estab- 
lishing angular dispersion and X-shaped ST spectra in 
the amplified modes [l2(. We expect that the same key 
action of the gain should play a role in our experiment 
as well. In fact, we detected appreciable angular dis- 
persion in the generated signal wave. Moreover, even if 
the temporal breakup in Fig. [21 (right) appears similar to 
what might be expected by a purely temporal MI, the 
numerical results in Fig. [31 indicate that the imbalance 
{i.e. Ep/F s <C 1) is not retrievable as long as spatial and 
temporal processes are considered independently. 

In order to provide an explanation of the observed con- 
nection between dimensionality, imbalance and depen- 
dence on n we focus attention on the dramatic increase 
of the down-conversion bandwidth close to degeneracy, 
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FIG. 2: Spatio-temporal intensity maps of the output pump- wave packet for OPA excitation and two different frequency-tuning 
conditions (top and bottom); left: experimental data; center: smoothed numerical data; right: high resolution numerical data. 
All data refer to the same input-pulse energy. 



which should cause a down-conversion rate much larger 
than that for up-conversion and thus the imbalance to ap- 
pear. Note that a correlation between bandwidths and 
imbalance must be expected as long as the signal and 
idler field are not correlated. In fact, in case of phase- 
conjugated fields, the expected drop in up-conversion ef- 
ficiency induced by the filtering action of the narrow up- 
conversion bandwidth no longer occurs, mutual cancela- 
tion of the opposite chirps ta king place without losses 
in the sum- frequency process [l3j. This consideration 
explains why the ID temporal model fails in retrieving 
the imbalance. In fact, a phase conjugation between sig- 
nal and idler is established by the down conversion pro- 
cess. However, in the real 3D space, this conjugation is 
quenched very rapidly owing to the fast phase accumu- 
lation caused by diffraction and angular dispersion [lif . 
In contrast, in the plane- wave approximation, it becomes 
unrealistically robust, due to the weaker effect of GVD. 

With the aim of verifying the generality of the pro- 
cess we performed a second experiment, where the spatial 
solitons are excited via SHG. The measurements, based 
on time-integrated detection only, were performed by fo- 
cusing down to a 45 /im spot (FWHM) a first-harmonic 



(FH) pulse (1.5 ps, 1055 nm) on the input face of a 
30 mm long LBO crystal, operated in non-critical, type-I 
phase matching. By changing the crystal temperature we 
tuned the phase mismatch parameter, Afc = 2fcpH — &sh, 
while keeping the frequency tuning fixed at degeneracy 
(f2 = 0.5). The results in terms of the ratio between the 
energy contents, ^fh/^sh, at threshold (both FH and 
SH beams trapped) are presented in Fig. 01 the Afc < 
region being limited close to phase matching by a rapid 
increase of the threshold with |Afc|. The few points, how- 
ever, are sufficient to highlight the large discrepancy oc- 
curring for Afc < between experiment and the CW- 
model, which predicts solitons with the opposite imbal- 
ance with respect to that measured. 

The fact that the CW model fails for negative, and not 
for positive, phase mismatch must be related to the dom- 
inant SH content expected in this parameter region, as 
long as spectral broadening does not occur (see Fig.EJ full 
line). Due to this, the dynamics that follows the SH gen- 
eration must resemble that of the OPA-excitation, where 
only the pump was launched at the input. To this end 
we mention that, for Afc < 0, a slight input-energy re- 
duction caused the quenching of the trapping only in the 
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FIG. 3: Frequency-tuning dependence of the pump-to-signal 
(peak) fiuence ratio, F p /F s , for OPA excitation. Squares: ex- 
perimental; Triangles: 3D+1 spatio-temporal model; Circles: 
2D+1 spatial model; Stars: 1D+1 temporal model. Fluences 
from the 3D+1 model are calculated by temporal integration. 
All data refer to the same input (peak) intensity. 
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FIG. 4: Phase-mismatch dependence of the ratio between the 
output FH and SH energies (experimental: balls), or intensi- 
ties (CW model: stars), integrated within the FWHM diam- 
eter. All data are at taken at threshold. 



SH wave; this means an identical red-soliton phenomenol- 
ogy as described for the OPA case. In comparing the two 
excitation schemes, we expect that for the SHG case clas- 
sical noise should also contribute to the process. As to 
the impact of the self-phase modulation, we mention that 
preliminary 3D+1 calculations performed in the absence 
of noise failed in retrieving the imbalance. A second pe- 
culiarity of the SHG excitation is that it leads to a true 
instability of the solitary regime, a spatial soliton being 
expected in the absence of MI. Again, since no theory has 
yet been developed for the MI of a bulk solitary beam, we 
must rely on the study of the MI of the FH+SH plane- 
and-monochromatic eigenmode for disclosing the mecha- 
nism of the ST interplay. Indeed, in this case too the MI 
has been shown to lead to the spontaneous formation of 



angular dispersion, its gain curve having an X-type pro- 
file in the ST frequency domain |l5j . 

In conclusion, we have shown that conventional, time- 
integrated, beam-profile measurements reveal the oc- 
currence of a x spatial soliton for which only the 
low-frequency signal and idler waves appear to be self- 
trapped. The implementation of a spatio-temporal, 3D 
mapping, technique has revealed that also the high- 
frequency pump is trapped, but with an unexpectedly 
weak intensity. Numerical calculations have highlighted 
the occurrence of a temporal break up of all the three 
fields on a lOfs scale. Moreover, they have shown that 
dispersion, diffraction and noise are all necessary ele- 
ments for describing the process. Relying upon recent 
studies on the spatio-temporal instability of the plane 
and monochromatic wave in OPA and SHG, we propose 
an interpretation of the process as the manifestation of 
a noise-seeded spatio-temporal modulational instability 
on the %( 2 ) spatial-soliton dynamics. 
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